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Abstract 

We obtain the induced Lorentz- and CPT- violating term in QED at finite temperature using 
imaginary-time formalism and dimensional regularization. Its form resembles a Cliern-Simons- 
like structure, but, unexpectedly, it does not depend on the temporal component of the fixed 
constant vector that is coupled to the axial current. Nevertheless Ward identities are respected 
and its coefficient vanishes at T = 0, consistently with previous computations with the same 
regularization procedure, and it is a non-trivial function of temperature. We argue that at finite 
T a Chern-Simons-like Lorentz- and CPT-violating term is generically present, the value of its 
coefficient being unambiguously determined up to a T— independent constant, related to the 
zero-temperature renormalization conditions. 



1 Introduction 

The phenomenological consequences of breaking Lorentz- and CPT-invariance have been actively 
explored in recent years as they could be measurable low-energy effects of quantum gravity or 
superstrings Q. In QED this issue was examined some time ago in Ref. ||], while lately CPT 
and Lorentz non-invariant extensions of the Standard Model were scrutinized in Ref. [Q] . As many 
breaking terms are allowed, most efforts have been focused on the possible constraints coming from 
experimental data H as well as from renormalizability requirements and anomaly cancellation. 
In this context, there arose a "theoretical" controversy on the possibility of generating, through 
radiative corrections, a Chern-Simons like term in the effective action of QED. There Lorentz and 
CTP symmetries can be in fact destroyed by considering a term of the form 

Cos = ^-e^A v F aP , (1.1) 

where k^ is a constant vector. This breaking term, suggested in Ref. ||], predicts birifrangence of 
the light in the vacuum and observations on distant galaxies put a very stringent bound on k^ (H . 
On the other hand the superstring inspired extensions of the Standard Model proposed in Ref. |4|] 
contain, in the fermionic sector, a Lorentz- and CPT-violating axial-vector coupling 

c b = (1.2) 

with a constant, prescribed four-vector that couples to the usual axial-current of QED. The 
interaction term Cb could generate, through radiative corrections, a non-vanishing value for k^ 
||. If this were the case, the strong bounds on k^ would translate into strong bounds for the 
non- invariant term (1.2). The aforementioned controversy arises from the fact the calculation is 



plagued by a dependence on the regularization adopted. While some papers claim that particular 
methods offer the correct result, others argue that the requirement of vector gauge invariance forces 



a vanishing induced term [||, ||. Recently this issue was also discussed in [[H]] in the heat kernel 
approach. A rather lucid discussion of the problem appeared in Ref. [11], where it was pointed 
out that the relevant form of the vectorial Ward identities may depend on how the vector 6„ is 



embedded into (or derived from) a more fundamental theory. As an example of that in [11] it was 
proposed an axion-like model to generate b^ as VEV of a dynamical field: there, a weaker form of 
the vectorial Ward identity governs the appearance of the interaction ( |1.1|) in the effective action, 
and, in particular, it is not strong enough to ensure the vanishing of its coefficient. In any case, 
when b^ is considered a strictly constant non-dynamical vector field, only vectorial Ward identities 
with vanishing axial momentum are relevant and they do not fix the actual value of the coefficient 
of the Chern-Simons term: it depends on the renormalization condition. 

Our interest is instead devoted to a different feature of the problem: the purpose of this letter 
is in fact to study the effect of a thermal bath on the structure of the Chern-Simons like term fll.lD • 
obtained by integrating out fermions coupled to the axial- vector ba- In particular our starting ob- 
servation is that regularization ambiguities cannot modify temperature dependence, since they are 
related to the ultraviolet behavior of the theory, that is temperature independent. Renormalization 
conditions are usually implemented at T = 0, where the parameters of the theory are defined, and 
consequently their temperature evolution is determined. Our one-loop computation may suffer, 
therefore, of the mentioned T = ambiguities while the functional form of the induced term and 
the temperature dependence of its coefficient are safe. We will use imaginary-time formalism and, 
for simplicity, dimensional regularization: in this scheme, where the vectorial Ward identities hold 
even at non-zero axial-vector momentum, a consistency check of our algebra is given by the van- 
ishing of the induced term ( |1 . 1| ) in the limit T —* 0. The fact that dimensional regularization in its 



standard form does not allow the appearance of the CS term at T = has been already pointed out 
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in || . On the contrary a CPT- and Lorentz- violating Chern-Simons action is generically present 
at T 7^ 0: while this fact may have some relevance for phenomenological application, potentially 
being active in the early universe, a serious question arises about the consistency of the effective 



theory. For time- like it was shown in [12] that the vacuum is unstable under pairs creation of 
tachyonic photon modes with finite vacuum decay rates and, recently, it was argued [JljJ that, in 
this case, unitarity itself may be in trouble (see also the original discussion in S). A more gen- 



eral analysis on the consistency of the theory at quantum level has been presented in [14], where 
both time-like and space-like cases appear to be problematic when microcausality and stability are 
examined. Rather surprisingly our computations show that the induced term does not depend on 
the temporal component of b^: we have not explored up to now the dynamical consequences of this 
fact in our finite-temperature context. Moreover invariance under small (i.e. not wrapping around 
the compactified imaginary time ) vectorial gauge transformations of the induced term is easily 
shown due to the use of dimensional regularization. 

2 The structure of one-loop self-energy at 

To begin with let us consider a modified QED action described by the Lagrangian density 

C = i>[i^-m-j 5 p-e4]ip. (2.3) 

As discussed in Ref. Q the linear contribution to the Chern-Simons term arises from the photon 
self-energy with one insertion of the axial-vector field, 

nf (p) = ib% luX ( P ) + Iu^x(-P)}, (2.4) 
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where I^x is given by the "triangle" -like graph with zero momentum axial vertex (we are working 
from now on directly in Euclidean space): 

T M- ■ 2 M 4fc ^Mff + wpTAT^ + imHM± t>) + im)} 

V * W Ze 7 (2vr) 4 + m 2 ) 2 ((A: + p) 2 + m 2 ) ' 1 j 

The CPT- and Lorentz- violating Chern-Simons action is extracted from Eq. ( |2.5| ) by isolating, from 

the odd-parity part, the tensorial structure linear in the external momentum and by performing 

the limit p 2 — ► in the scalar integral multiplying it. In Ref. B the explicit evaluation at T = of 

(p) in the dimensional regularization was presented. In particular it was noticed that the only 

algebraic properties of 75 used in the computation were: a) the trace of 75 with an odd number 

of Dirac matrices vanishes and b) the trace of 75 with an even number of Dirac matrices can be 

reduced using the Clifford algebra to the quantity Tr[7^7i/7a7 j 375]- Consistency also requires that 

Tr[ 7/ ,7,7 5 ] = TV[ 75 ] = 0. 

In this zero-temperature case, the linear p^ dependence is easily extracted and the result can 

be presented as 

e 2 r i 

" ^ fe VTr [ 7 ^7A7/375] [Fi{p 2 /m 2 ) + F 2 {p 2 /m 2 )\ , (2.6) 

where the explicit form of J-\(p 2 /m 2 ) and J-2(p 2 /m 2 ) is given in Q: evaluating J-i(p 2 /m 2 ) and 
J-2(p 2 /m 2 ) in D dimensions, taking the limit D = 4 and expanding in p 2 , it results that 

T^jm 2 ) + F2(p 2 /m 2 ) ~ (2.7) 

showing the absence of a Chern-Simons contribution to the effective action. We agree with this 

computation but we want to recall a couple of remarks in order to better appreciate the finite T 

effects. First of all one can check that the cancellation of the leading order (constant) contribution 

to J-\{p 2 1 m 2 ) + J-2(p 2 /m 2 ) comes from a delicate balance between a "classical" term (proportional 

to m 2 in Eq.( |2.5[ )) and an "anomalous" quantum term, deriving from the potential divergences^. 
1 This mechanism of cancellation of a quantum contribution against a classical one is reminiscent of the analogous 



As mentioned in Ref.[ll], the massless case escapes this mechanism, therefore dimensional regular- 
ization gives there a non- vanishing result. This fact is related to the loss of analyticity || of the 
"triangle" -like diagram at m = in the limit of vanishing axial- vector momentum. In the finite 
temperature case, where the analyticity properties in the external momenta are usually weaker, 
this observation suggests the concrete possibility that a non-zero result could appear even in the 
massive case. Secondly, we stress that the covariance of the momentum integration immediately 
selects the Chern-Simons tensorial structure and the dependence on p 2 of the coefficient function 
at zero-temperature. This is no longer true at finite temperature, as we shall see in a while, due to 
explicit presence of the Matsubara frequencies: in particular the double limit po — ► and p% — > 
has to be performed very carefully. 

Let us assume, from now on that the system is in thermal equilibrium with a temperature 
T = /3 _1 , (3 being interpreted as the radius of the compactified euclidean time. In this case 
we may use the Matsubara formalism that consists simply in taking ko = (n + (anti-periodic 
boundary conditions for fermions requires semi-integers frequencies) and replacing ^ / dko = 4 Y^n- 
The remaining J d 3 k integral is of course continued to D-spatial dimensions. The trace can be still 
performed in full generality and simple algebraic manipulations in the loop momenta (not involving 
shifts or symmetry properties) allows us to write fl2.8|) as follows: 

Va(p) = 

_4-^V f d ° k ( 2fc M^Ap CT K^-2^e MpCT p^)-2(A:-p)e^ A X+p 2 e^ Ap (^-p^ 
l (3^J{2n) D {k 2 + m 2 ) 2 {(k+p) 2 +m 2 ) 

f3^J (2tt)D (k 2 + m 2 ) 2 ■ ^ 



phenomena in 3D in the case of parity anomaly. This perfect balance between the two contributions is peculiar of 
the standard dimensional regularization. In other scheme this exact cancellation does not occur, leaving us with a 
non vanishing CS term whose coefficient is however temperature independent. 
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The following step in the computation is to extract the tensorial structure, leaving us with the 
evaluation of scalar integrals. The second term in Eq, (|2.8| ) is easily tamed, both the integral over 
the spatial components and the series over ko are anti-symmetric in exchanging k — > —k (we can 
find a region around D = 3 where everything is convergent). It remains therefore 



e 2 



- 4i— e tlu x P P P Io, (2.J 



where we have 



that exhibits the Chern-Simons like structure. Let us discuss now the first contribution. We 
introduce Feynman parameters in order to perform the integral. To implement translations only 
on the space components of the loop momentum we decompose k^ as follows 

hp = hp + ko8 0fJl . (2.11) 

Shifting k — * k — xp (where p^ is defined as above) in Eq.(^j^) and using the covariance under 
spacial rotations, which allows us to conclude that 

A , k 2 

we arrive to the form 

e 2 

I/iux(p) = _4d -g [ e iMu\oh + e^x P P p (^h + h) 
+ (p p e uXp oP p ~ P v £ tl \ P oP p ~ P 2 £pv\o)h 

+ (^e l/Xp oP p ~ <We M A P oP p ~ Poe^AoKJ , (2-12) 



where 



1 , o/, ^ f d D k [p 2 {l-2x)(k + xp ) - %k 2 p + 2p {k + xp ) 2 



J q dx2(l-x)J2 J 



(2tt) D [^2 _|_ ^ Q _|_ X p Q )2 _|_ x (l _ x ^p2 _|_ m 2]3 



dx2(l-x)T I [-^o + ^-^-M 1 -^] 

o >Z ^J (2tt) d + {k Q + xpo) 2 + x{1 - x)p 2 + m?Y' 



d D k 



2x(k + xpo) 



(2it) d [k 2 + (ko + xpo) 2 + x(l - x)p 2 + m 2 ] 



dx2(l-x)J2 



d D k 



2 t2\ 



(2(k + x Po ) 2 - ±k 



(2.13) 



(2tt) d p + (k + xp ) 2 + x(l - x)p 2 + m 2 f 
Working at finite temperature, more structures have been generated, a fact that is not unexpected 
due to the explicit breaking of four-dimensional covariance. The important point is that, nev- 
ertheless, the tensors must be transverse respect p^ and p u , since the vectorial Ward identity are 
unaffected by the presence of the temperature. By inspection we see that the only potential trouble 
comes from I\ (being e^ v \Q not transverse!). Luckily we can show that I\ is exactly zero. To this 
purpose, it is useful to rewrite I\ as follows 



2p k 2 /D 



k + xp d 
2 dx 



1 



[k 2 + (£;o + xpo) 2 + x(l — x)p 2 + m 2 ] 2 

(2.14) 



[k 2 + (ko + xpo) 2 + x(l — x)p 2 + m 2 ] 3 
and integrating by part with respect to x (the boundary terms are zero) we get, after having 
performed the D— dimensional integral, 

(k + xpo) 



(4tt) 2 Jo [(ko + xpo) 2 + x(l — x)p 2 + m 2 

We can use now the explicit form of the Matsubara frequencies and the fact that po is discrete 



,212-S 



(2.15) 



(po = relabeling the sum in Eq.( 2.15 ) as n —>■ —n — 1 — I and making the change of variables 



y = 1 — x, one easily obtains that 



ii 



-h. 



We stress that we did our computations in D dimensions, where everything is convergent and no 
limit on p has been performed. The same arguments applies to ^3, and we remain, therefore, 
with two independent tensorial structures: we need to evaluate Iq + 2I2 and I4. Before entering the 



computations we remark that the emergence of a new, transverse tensorial structure was overlooked 
in Ref. flGfl , where the coefficient of the Chern-Simons term at finite T was obtained simply 
evaluating the scalar integral, relevant at T = 0, by introducing Matsubara frequencies for pq. As 
we will see in the next section, our result disagrees with that. 



3 The CPT- and Lorentz- violating term at T ^ 



Let us evaluate the coefficients of the two independent structures, in the small momentum limit: 
we remark that at finite temperature this procedure is rather delicate, due to the fact that, in 
general, the limits po — ► and p 2 — > do not commute [17]. Here we shall take first pq — > and 
send p 2 — ► 0. The sum Iq + 2I2 becomes, when the D-dimensional integrals have been computed: 



Jo + 2J 2 



r(2-S] ,! 



D , dr2(l-x)£ 
(4vr)- Jo 



1 



D 



1 



[k 2 + 



in 



212- 



+ 



4 2-D/2 



1 



J^ r(3) [fcg + m 2 ] 2 - 



(3.16) 



The result is zero identically in arbitrary dimension: one can check that when the dependence on p 2 
is retained the corrections are regular, and of order p 2 (and of course the coefficient depends on T). 
We see that the the zero-temperature tensorial structure still has a vanishing coefficient in the small 
momentum limit when The only possible contribution to Lorentz- and CPT-violation could 

therefore arise from I4, i.e. from the non-covariant structure. The relevant term to be calculated is 



Ji 



2 l2\ 



d D k (2feg - %k 
(2vr) D p + k 2 + m 2 } 



(3.17) 



The D-dimensional integration leads to 



T(2 



D 

2 . 



(4vr 



dx(l — x) 



(3-D) 



[kg + m 2 ] 2 2 



+ (D-4) 



m 



[k 2 + m 2 ] 3 2 



(3.18) 
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At this point we need an explicit representation for the sum over the Matsubara frequencies: we 



use the following result [18], valid when 1/2 < A < 1 



E[(n + ^ + a T A = ^;^ 2 ' + 4sin( T A) 



dz 



-Re 



| a | (z 2 - a 2 ) x \exp 27r(z + ib) — 1 



. (3.19) 



We cannot apply directly this formula to our case: at the end we want to take the D = 3 limit, 
and it is clear that, evaluating the second contribution to I4 in Eq,( |3.18 ), the integral in Eq.( 3.19 ) 
do not converge there in the limit D = 3. It is not difficult anyway to perform the analytical 
continuation in Eq,( 3.19| ) using the relation 



dz 



\ a \ (z 2 - a 2 ) x 



Re 



1 3-2A f°° 



exp 2ir(z + ib) — 1/ 2a 2 1 — A J\ a \ 



dz 



(z 2 - a 2 ) x - 1 



Re 



1 



dz 



Aa 2 (2 - A)(l - A) J\ a \ {z 2 - a 2 ) x ~ 2 dz 2 



Re 



1 



exp 2ir(z + ib) — 1 



exp 2ir(z + ib) — 1, 

(3.20) 



Eq.(3.18) can now be explicitly evaluated at D = 3: we see that the potential singularity at D = 3, 
coming from the first contribution, cancels (notice the factor D — 3 in front), the finite residue 
(that would be temperature independent) cancels with an analogous term coming from the second 
contribution, leaving us with the final result 



/■OO 

I 4 = 2/3 / dz(z 2 - £ 
J\£\ 



2 t 2\i tanh(-7T2:) 



cosn (ttz) 



(3.21) 



where we have defined ^ = The behavior of F(£) is displayed in the plot below. Eq.( 3.2l| ) is 




Figure 1: plot of F(£) 
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the main result of our paper: it shows that for (3 ^ oo (T ^ 0) a CPT- and Lorentz-violating term 
appears. In momentum space it can be written as 

nf(p) = ±e 2 F(Ob x (5 0fM e uXp oP p - 5 u^x P oP p - Po^uxo) + 0(p 2 ). (3.22) 

Several comments are now in order: first of all when T = F(£) vanishes, recovering therefore the 
fact that, using dimensional regularization, no CPT- and Lorentz-violating Chern-Simons like term 
is present in the effective action. The opposite limit (T — > oo) is otherwise finite (F(0) = l/2n 2 ): 



in Ref. [16] a similar behavior was found for the temperature evolution of the coefficient of the 
Chern-Simons like term but there the T = boundary condition was taken so that at T = oo 
the symmetries were restored. Moreover, at variance with our result, the Chern-Simons term 
there was implicitly assumed to be related to the covariant tensorial structure, fact that from our 
computation turns out to be incorrect. The second point is that our induced action does not depend 
on the temporal component of b^; in refs. [12, 13, 14 1 it was discussed the consistency of the 
theory at quantum level, when the Chern-Simons like action is present. It would be interesting to 
address the problem of stability in the finite-temperature situation considering our induced term. 
Another observation is related to the dependence on £: we see that the limit (3 — > is the same as 
m —* (since the only dependence on the mass and on the temperature appears through £). This 
suggests that the presence at finite temperature of a non-zero CS-like term is related to the loss 
of analyticity in external momenta, bypassing therefore the argument of Coleman and Glashow || 
against it (analyticity was also assumed in Ref. pi) . It is interesting to write down the induced term 
in configuration space 

S«cs» = 4ie 2 F(C) J d A x bi[A^ k F jk - 2e i ^ k A J F 0k ] (3.23) 
gauge invariance is achieved via Bianchi identity up to a total derivative. 
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The asymmetrical behavior played in the above action by the "spatial" and "temporal" com- 
ponent of bfj, (the latter being completely absent) might seem strange. A source of this asymmetry 
can be surely traced back to the presence of a thermal bath which, selecting a specified frame, 
provides an additional Lorentz violation. However, below, with the help of the analogous problem 
in two dimensions, we would like to suggest that the origin of this term at finite temperature may 
be related to a deeper geometrical reason. In D=2, in fact, its structure and its coefficient can be 
easily understood through the interplay of the global part of the Quillen anomaly, which controls 
the obstruction to the chiral splitting, and the presence of a non trivial cycle (finite temperature 
in the above language). 

There, the relevant Green function is the 1-point function with one 6^-insertion (we study the 
massless case for simplicity). Using dimensional regularization we see that the zero-temperature 
computation gives a vanishing result because / d D k^^ = (we remark that having constant 
implies that the external momentum has to be null) . But turning on the temperature the situation 
changes drastically. The relevant integral is: 

d D k k a k u 



D* = - ie b x Tr[ 7 , 7a7m ]W J 4^D^> ( 3 - 24 ) 



that is equivalent, after using Dirac algebra to 



2eV 



i r d D % i 1 ^ f d D k fcg 



l 3 (2tt)D [k 2 + k 2 } wu " w ^i (2vr)^ p + k{ 

+9 J. V f d °" k 
+ Xa 0^J (2^) D \kl + k 



212 



(3.25) 



'P^J (27r)^[^ + A;2]2_ 

It is not difficult to see that the D— dimensional integration and the analytical continuation of the 
sum gives a finite result for the b± component (due to a cancellation between a pole and a zero as 
D = 1) leaving us with 

n£ = — eb x e Xa 5 ^ (3.26) 
12 



while one can easily show that 60 component has zero coefficient after performing the ^-integration. 
We remark that here no external momentum limit has been done, therefore the result is exact. The 
above term has a natural interpretation in configuration space as 



i J (fx^Aoix), (3.27) 



that is the analogue of D = 4 (the "complete" CS-like term would be here e^b^Ay). The term 
appearing eq. (|3.27| ) is nothing but a remnant of the holomorphic anomaly on the torus [19|. In 



fact, when non-trivial cycles are present in the base-space, the effective action acquires a subtle 
dependence on the harmonic part of the gauge potentials. In particular when both vector and axial 
gauge fields are coupled, requiring gauge invariance implies that an anomalous phase has to be 
present in order to cure the transformation properties of the modulus of the Dirac determinant. This 
phase can be derived from general algebraic-geometrical arguments [^] , being related to the Quillen 
anomaly, or by an explicit £— function computation [gCJ] of the relevant functional determinants. 
The asymmetric character of the phase has to be ascribed to the anomalous modular transformation 
properties of chiral partition functions. We can now understand the appearance of this term at 
finite temperature: b^ is basically an harmonic one-form axially coupled and therefore being able to 
interact with the harmonic component of Aq. The complete anomalous phase requires a quadratic 
part in 6 M , that can be easily recovered by computing the Feynman graph with two 6 M insertions. 
This discussion may suggest that the four-dimensional term could have an appealing mathematical 
interpretation. 
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